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Abstract By using realistic models for elastic hadron scattering we demonstrate that at 
; present accelerator energies the s-channel unitarity bound is safe, not to be reached until 

10 5 GeV, while the black disc limit is saturated around 1 TeV. It will be followed by a larger 
■^j- ■ transparency of the scattered particles near the center. 



1 Introduction 



o 

On 
On 

^3 
Oh! 

Our decision to write this paper was motivated partly by recent claims that in high energy 
hadron scattering the black disc limit has been reached and the violation of the s-channel 
unitarity in some models is just around the corner. While the first statement is true and has 
interesting physical consequences, the second one is wrong for any realistic model fitting the 
existing data on proton and antiproton scattering up to highest accelerator energies. 

To start with, let us remind the general definitions and notations. 
Unitarity in the impact parameter (b) representation reads 

%mh(s,b) = \h(s,b)\ 2 + G m (s,b), (1) 

where h(s, b) is the elastic scattering amplitude at yfs center of mass energy (with 3m/i(s, b) 
usually called the profile function, representing the hadron opacity) and Gi n (s, b), called the 
inelastic overlap function, is the sum over all inelastic channel contributions. Integrated 
over b, (1) reduces to a simple relation between the total, elastic and inelastic cross sections 
o-tot(s) = a e i(s) + a in (s). 

Equation (1) imposes an absolute limit 

< \h(s,b)\ 2 < Smh(s,b) < 1, (2) 
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while the so-called "black disc" limit a e i(s) = o-, in (s) = |cr tot (s) or 

Sm/i(s, b) = 1/2 (3) 

is a particular realization of the optical model, namely it corresponds to the maximal ab- 
sorption within the eikonal unitarization, when the scattering amplitude is approximated 

as 

fc(s,&) = ^(l-exp[ia;(s,b)]), (4) 

with a purely imaginary eikonal u(s, b) = iQ(s, b). 

Eikonal unitarization corresponds to a particular solution of the unitarity equation 



Ks,b) = \ 



1 ± J l - AG in (s, b) 



(5) 



the one with minus sign. 

The alternative solution, that with plus sign is known [ffl, [| and realized within the 
so-called [/-matrix approach ||, |J] where the unitarized amplitude is 

where now U is the input "Born term", the analogue of the eikonal us in (4). 

In the [/-matrix approach, the scattering amplitude h(s, b) may exceed the black disc limit 
as the energy increases. The transition from a (central) "black disc" to a (peripheral) "black 
ring", surrounding a gray disc, for the inelastic overlap function in the impact parameter 
space corresponds to the transition from shadowing to antishadowing We shall present 
a particular realization of this regime. 

The impact parameter amplitude may be calculated either directly from the data, as it 
was done e.g. in || |6| (where, however, the real part of the amplitude was neglected) or 
by using a particular model that fits the data sufficiently well. There are several models 
appropriate for this purpose. In a classical paper [|7j on this subject, from the behaviour 
of Gi n (s,b), the proton is characterized as getting "BEL" (Blacker, Edgier and Larger). As 
anticipated in the title of our paper, the proton, after having reached its maximal darkness 
around the Tevatron energy region, may get less opaque beyond. 

Actually, the construction of any scattering amplitude rests on two premises : the choice 
of the input, or "Born term" and the relevant unitarization procedure (eikonal or [/-matrix 
in our case). Within the present accelerator energy region there are several models that 
fit the data reasonably well. Compatible within the region of the present experiments, they 
differ significantly when extrapolated to higher energies. We shall consider two representative 
examples, namely the Donnachie-Landshoff (D-L) model ||, ^| and the dipole Pomeron (DP) 
model 1,0) 

In Sec. 2 we present the necessary details about the two realistic models (D-L and DP), 
then, focusing on the DP model, we investigate in Sec. 3 the unitarity properties at the 
" Born level" and in Sec. 4 we study the optical properties (transparency) after unitarization; 
a comparaison with the D-L model is given in the Appendix. 



4 We follow traditional terminology, although the word "matrix" in this context is misleading, since U, 
similar to the eikonal, is a single function rather than a matrix. 
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2 The "Born term" 

The Donnachie-Landshoff (D-L) model || is popular for its simplicity. Essentially, it means 
the following four-parametric empirical fit to all total hadronic cross sections 

a tot = X s 5 + Y s s % (7) 

where two of the parameters, namely 5 = «p(0) -1« 0.08, and S r (< 0) are universal. While 
the violation of the Froissart-Martin (F-M) bound, 

°tat{s) < C (£ns) 2 C = 60 mb , (8) 

inherent in that model, is rather an aesthetic than a practical defect (because of the remote- 
ness of the energy where eventually it will overshoot the F-M limit), other deficiencies of the 
D-L model (or any other model based on a supercritical Pomeron) are sometimes criticized in 
the literature, but so far nobody was able to suggest anything significantly better instead. A 
particular attractive feature of the D-L Pomeron, made of a single term, is its factorizability, 
although this may be too crude an approximation to reality. 

The t dependence in the Donnachie-Landshoff model is usually chosen || in the form 
close to the dipole formfactor. For the present purposes a simple exponential residue in the 
Pomeron amplitude will do as well, with the signature included 

A(s,t) = -N (-Z-) e Bt , (9) 
V s dl J 

where a(t) — a(0) + oi t is the Pomeron trajectory and N is a dimensionless normalization 
factor related to the total cross section at s — s<$ by the optical theorem 

N= . Sd [ M cr tot (s = s dl ) . (10) 
Air sin |a(0) 

According to the original fits g, g|: s di = 1 Gey2 > «(°) = 1 - 08 > «' = °- 25 GeV~ 2 , and 
X = 21.70 mb (see (7)) resulting in iV = 4?r sin ^a(o)/2 = 4.44. By identifying 

da (s,t) _ da{s,t = 0) cBexv(s) t 
dt dt 

and choosing the CDF or E410 result for the slope B exp at the Tevatron energy, we obtain 
B = \B exp (s) - a'in^ = 4.75 GeV" 2 . 

In the dipole Pomeron (DP) model |4j], factorizable at asymptotically high energies, 
logarithmically rising cross sections are produced at a unit Pomeron intercept only and thus 
DP does not conflict with the F-M bound. While data on total cross section are compatible 
with a logarithmic rise (DP with unit intercept) the ratio o e ijo tot is found (see [0 for 
details) for 5 = to be a monotonically decreasing function of the energy for any physical 
value of the parameters. The experimentally observed rise of this ratio can be achieved 
only for 5 > and thus requires the introduction of a "supercritical" Pomeron, a(0) > 1. 
As a result, the rise of the total cross sections is driven and shared by the dipole and the 
"supercritical" intercept. The parameter 5 = a(0) — 1 in the DP model is nearly half that 
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Table 1: Parameters of the Dipole Pomeron found in 



a 


b p 


Qt(0) 


a'(GeV- 2 ) 


e 


s (GeV 2 ) 


355.6 


10.76 


1.0356 


0.377 


0.0109 


100.0 



of the D-L model, making it safer from the point of view of the unitarity bounds. Generally 
speaking, the closer the input to the unitarized output, the better the convergence of the 
unitarization procedure. 

Let us remind that apart from the "conservative" F-M bound, any model should satisfy 
also s-channel unitarity. We demonstrate below that both the D-L and DP model are well 
below this limit and will remain so within the forseable future. (Let us remind that the 
D-L and the DP model are close numerically, although they are different conceptually and 
consequently they extrapolations to superhigh energies will differ as well.) 

The elastic scattering amplitude corresponding to the exchange of a dipole Pomeron reads 



A(s,t)= A [e-W2 G(a )( s/so ) 

= e- ina / 2 {s/s ) a [G'{a) + (L — m/2)G(a)} 



(12) 



g 

where L = In — , a = a(t) is the Pomeron trajectory; in this paper, for simplicity we use a 
s 

linear trajectory a(t) = a(0) + a't. 

By identifying G'{a) = -ae b ^ a ~ l \ (12) can be rewritten in the following "geometrical" 
form 

A( s , t) = i^- \r 2 x {s) e^MM')- 1 ! - e r 2 2 {s) e r ^ a ^} , (13) 

where 

rl(s)=bp + L-i^ , r \{s) = L-i X -. (14) 

The model contains the following adjustable parameters: a, b p , a(0), a', e and so- 

In Table 1 we quote the numerical values of the parameters of the dipole Pomeron fitted 
T0|| to the data on proton-proton and proton-antiproton elastic scattering : 



in 



a tot (s) = -ZmA(s, 0) , p(s) = I^T^ i 4 < v^GeV) < 1800 (15) 
as well as the differential cross-section 

da(s,t) = tt_^ s ^2 . 23.5 < v^(GeV) < 630 ; < |t|(GeV 2 ) < 6 . (16) 

(Xb S 

In that fit, apart from the Pomeron, the Odderon and two subleading trajectories u and / 
were also included. Here, for simplicity and clarity we consider only the dominant term at 
high energy due to the Pomeron exchange with the parameters fitted in [|TJJ. The extent to 
which this Pomeron is a good approximation in the TeV region is discussed in details in [[12] . 



The quality of this fit is illustrated and discussed in [10]. With such a simple model and 
small number of parameters, better fits are hardly to be expected. 
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We use the above set of parameters to calculate the impact parameter amplitude, and 
to scrutinize in Sec. 3 the unitarity properties of this "Born level" amplitude. In Sec. 4 we 
introduce a unitarization procedure, necessary at higher energies and discuss the relevant 
physical consequences. 

To summarize, the DP model with a unit intercept is selfconsistent in the sense that its 
functional (logarithmic) form is stable with respect to unitarization. Moreover, the presence 
of the second term, proportional to e in (13) has the meaning of absorptions and it is essential 
for the dip mechanism. It can be viewed also as one more unitarity feature of the model. 
In the limit of very high energies, when L ^> b p the two (squared) radii R\ = a'rf become 
equal and real and the model obeys exact geometrical scaling as well as factorization (see 
next section). Alternatively, it corresponds to the case of no absorptions (e = 0). 
However attractive, the case of a unit intercept (5 = 0) is only an approximation to the more 
realistic model, requiring 5 > to meet the observed rise of the ratio a e i/&tot- For such a 
"supercritical" Pomeron unitarization becomes inevitable. 



3 Impact parameter representation, the black disc limit 
and unitarity 

The elastic amplitude in the impact parameter representation in our normalization is 

1 r 00 

h(s, b) = — / dq qJ (bq)A(s, -q 2 ) , q = ^ . (17) 

AS JO 

The impact parameter representation for linear trajectories (] is calculable explicitly for 
the DP model (13) 

h(s,b) = i g [e r ' s e~ h2 / m " - e e r * 5 e^ 4 ^] , (18) 

where 

R 2 = a'r 2 (i = l,2) ; g = . (19) 
Asymptotically (i.e. when L ^> b p , i.e. y/s 3> 2. TeV, with the parameters of Table 1), 

h(s,b) g(s) (l-e)e-3& , (20) 

where 

R 2 = a'L ; 9(s) = 9o ■ (21) 

To illustrate the s-channel unitarity, we display in Fig. 1 a family of curves showing the 
imaginary part of the amplitude in the impact parameter-representation at various energies; 
also shown is the calculated (from (1)) inelastic overlap function. 



3 Other cases were treated e.g. in 
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Table 2: Central opacity of the nucleon 5smh(s, 0) calculated at ISR, SPS, Tevatron energies 
compared with experiment. 





53 GeV 


546 GeV 


1800 GeV 


exp 
th 


0.36@ 
0.36 


0.420 ±0.004 [I3| 
0.424 


0.492 ± 0.008 [U\ 
0.461 




12 3 12 3 

b (fm) 

Fig. 1. Calculated "Born level" ^smh(s, b) and Gi n (s,b) plotted versus the modulus of 
the impact parameter b for some characteristic energies yfs as indicated (solid curve for 
the LHC energy). The top of the scale on the left is the unitarity limit and the value 1/2 
corresponds to the black disc limit. The calculations are performed for the dipole Pomeron 
model; similar results are obtained for the D-L model, see the text. 

Our confidence in the extrapolation of ^smh(s, b) to the highest energies rests partly on 
the good agreement of our (non fitted) results with the experimental analysis of the central 
opacity of the nucleon (see Table 2). 

It is important to note that the unitarity bound 1 for 3m/i(s, b) will not be reached at 
the LHC energy, while the black disc limit 1/2 will be slightly exceeded, the central opacity 
of the nucleon being Qmh(s, 0) = 0.54. 

The black disc limit is reached at y/s ~ 2 TeV, where the overlap function reaches 
its maximum -r. This energy corresponds to the appearance of the antishadow mode in 
agreement with the general considerations in [|TJ. Notice that while $smh(s,b) remains 
central all the way, Gi n (s, b) is getting more peripheral as the energy increases starting from 
the Tevatron. For example at \/s = 14 TeV, the central region of the antishadowing mode 
below b ~ 0.4 fm is discernible from the peripheral region of shadowing scattering beyond 
b ~ 0.4 fm, where Gi n (s,b) = |. In terms of J7|, the proton will tend to become more 
transparent at the center ("gray", in the sense of becoming a gray object surrounded by a 
black ring), i.e. it is expected to become "GEL" instead of "BEL". 
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The s channel unitarity limit will not be endangered until extremely high energies (10 5 
for the D-L model and 10 6 GeV for the DP), safe for any credible experiment. It is interesting 
to compare these limit with the limitations imposed by the Froissart-Martin bound: actually 
the Pomeron amplitude saturates the F-M bound at 10 27 GeV. As expected, the F-M bound 
is even more conservative than that following from s-channel unitarity. 

The D-L and DL models are confronted in the Appendix. 



4 Unitarization 

Now, we consider the unitarized amplitude according to the "[/-matrix" prescription || |J 

with the "Born term" h(s,b) defined in the previous section in (13)-(14). 

Fig. 2 shows the behavior of the unitarized impact parameter amplitude H(s, b) and the 
corresponding inelastic overlap function at various energies. By comparing it with similar 
curves (Fig. 1) obtained at the "Born level" we see that unitarization lowers significantly 
both the elastic and inelastic impact parameter amplitudes. 




b (fm) 

Fig. 2. Same as in Fig. 1, for the unitarized amplitude H (s, b) and the overlap function, 
calculated without refitting the parameters used at the Born level. 

An unescapable consequence of the unitarization is that, when calculating the observ- 
ables, one should also replace the Born amplitude A(s,t) by a unitarized amplitude A(s,t) 
defined as the inverse Fourier-Bessel transform of h(s, b) 

A{s,t) = 2s dbbJ (by/^i)H{s,b) . (23) 
jo 
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Thus, the above picture may change since the parameters of the model should in principle 
be refitted under the unitarization procedure (this effect of changing the parameters was 
clearly demonstrated e.g. in [Hp . 

Actually, searching for a new fit of the parameters using an unitarization procedure is 
time consuming and unnecessary for the present discussion because the behavior of the 
amplitude and of the overlap function in the impact parameter representation obtained at 
the Born level will be almost restored after unitarization. We checked that the parameters 
of the complete model (with the secondary Reggeons and Odderon added in the fit) after 
unitarization may be rearranged so as to reproduce well the data and give roughly the same 
extrapolated properties as at the Born level. 

While the unitarity limit now is secured automatically (remind that Sm/i(s, 0) is well 
below that limit even at the "Born level" in the TeV region) the behaviour of the elastic 
impact parameter amplitude after it has reached the black disc limit corresponds (see [p]]) 
to the transition from shadowing to antishadowing. In other words, the proton (antipro- 
ton) after having reached its maximal blackness around 2 TeV, will become gradually more 
transparent with increasing energies at its center. 



5 Conclusions 

While the results of our analyzis in the impact parameter representation are in agreement 
with the earlier observations that Qrn/i(s, b) remains central and Gj n (s, b) becomes peripheral 
as the energy exceeds 2 TeV (see Fig. 1). There is a substantial difference with the known 
"BEL-picture" J7[, according to which with increasing energy the proton becomes Blacker, 
Edgier and Larger. 

We confirm that getting edgier and larger, the proton, after reaching its maximal black- 
ness, will tend to be more transparent or ("GEL", i.e. a gray disc surrounded by a black 
ring) when the energy exceeds that of the Tevatron. This transition from shadowing to new 
antishadowing scattering mode is expected to occur at the LHC. 

To conclude, we stress once more that both the data and relevant models at present 
energies are well below the s-channel unitarity limit. In our opinion, deviations due to the 
diversity of realistic models may result in discrepancies concerning 9 ; m/i(s, 0) of at most 10%, 
while its value at 2 TeV is still half that of the unitarity limit, so there is no reason to worry 
about it! Opposite statements may result from confusion with normalization. Therefore, 
model amplitudes at the "Born level" may still be quite interesting and efficient in analyzing 
the data at present accelerator energies and giving some predictions beyond. The question, 
which model is closer to reality and meets better the requirements of the "fundamental 
theory" remains of course topical. 

Extrapolations and predictions to the energies of the future accelerators (see e.g. ]T2| ) 
are both useful and exciting since they will be checked in the not-so-far future at LHC and 
other machines. The fate of the "black disc limit" is one among these. 

Acknowledgements We thank S.M. Troshin for a useful correspondance. 
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Table 3: Maximum values of the amplitude and overlap function at the Born level and after 
[/-matrix unitarization calculated at 14 TeV for the DP and D-L models without refitting 
the parameters 





Qmh(s, 0) 


G in (s, 0) 


SmH(s, 0) 


G in (s,0) 


DP 


0.535 


0.247 


0.349 


0.227 


D-L 


0.539 


0.246 


0.351 


0.227 



APPENDIX 

Comparaison between the DP and D-L models 

The D-L amplitude in the impact parameter representation at the Born 
level, calculated from (9) and (17) is 



h(s : b) = --(-i-) — — B'(s) = B + a f (£n--z-). (24) 
2s V s d iJ 2B'(s) s d i 2 

As already noted, the s channel unitarity limit both for the DP and the D- 
L model will not be endangered until extremely high energies (10 5 GeV for 
the DP and 10 6 GeV for the DP model, the order-of-magnitude differences 
coming from the smaller intercept in the DP model), while the F-M bound 
is saturated at 10 27 GeV (for more details see 

Table 3 presents a selection of results concerning the DP and D-L mod- 
els for the Pomeron in the impact parameter representation of the elastic 
amplitude and inelastic overlap function, calculated at b = at the LHC 
energy. 

We conclude that the two models give similar results; all conclusions on 
unitarity and black disc limits for DP model hold for D-L model as well (the 
curves in Figs. 1,2 would be indistinguishible by eye). 

Note that both models are supercritical, with asymptotic s s type behavior 
of the total cross sections. They are known to give fits which cannot be 
discriminated by present data from an asymptotic £n 2 s type behaviour. This 
is another argument to neglect unitarization effects. 
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